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ABSTRACT 

It is shown that a generalization of the use of “Vaisala orbits”, briefly mentioned 
at the Asteroids II conference, can be very conveniently accomplished by means of an 
inversion of the “GEM” form of the Gauss method. The procedure can also be applied to 
Apollo objects and to indeterminate cases of normal three-observation orbit computation, 
and there is also a simple extension to situations involving four or more observations. 


A popular procedure for planning and then identifying the third night’s observations of 
an asteroid is to utilize sets of orbits fitted to positions on the first and second nights on the 
assumption that the object is then at perihelion — or aphelion, a circular solution being the 
transition between the two sets. In the proceedings of the Asteroids II conference Bowell et al 
(1989) briefly considered that this procedure (Vaisala 1939) might be generalized to delineate 
the precise region of the sky in which a main-belt asteroid must lie. In a typical application 
the third night will be in the dark of the moon following that of the first and second nights, 
which — to avoid incorrect linkage — should themselves be separated by no more than about 
five days; and the two lunations will tend to surround or adjoin the asteroid’s being at 
opposition. Here I shall discuss in more detail how the “generalized Vaisala procedure” can 
conveniently be carried out in practice. 

If the asteroid had actually been identified on the third nijsjht, it would be reasonable to 
attempt a three-observation orbit computation from one position on each night. This could 
be accomplished by the Gauss (1809) method, preferably in the rigorous “GEM” (Gauss- 
Encke- Merton) form I described (Marsden 1985) in response to Taff’s (1984) misguided claim 
that the Gauss method is not mathematically valid. 

In the GEM form the orbit computation is made in Cunningham’s (1946) coordinate 
system, where the xy plane passes through the observations at times and £ 3 , the x axis 
being directed toward that at £ 3 . For reasons that will shortly become apparent, I shall 
consider 1 1 to refer to the first observation and £3 to the second observation. The time t 0 will 
then refer to the hypothetical third observation. The Cunningham coordinate system has 
the great advantage that the usual vector equation of the Gauss method separates into 

Po = (Zo - C\Z\ - c 3 Z 3 )/n 0 

Pi = (po™o - Yo + C 1 Y 1 + c z Y z )/c l m 1 (1) 

Pz = (po^o - £\Pih - Xo + c\Xi + c 3 X 3 )/c 3 , 

where the ( ) and (X t , l^,Zi) represent the components of the unit vectors from the 
observer to the asteroid and of the vectors from the observer to the sun at the times The 
scalar distances pi from the observer to the asteroid at the U are initially unknown, as are 
the scalars C\ and c 3 , the usual ratios of the areas of the triangles that are associated with 
the Gauss method. Eqs. ( 1 ^) can simply be evaluated in turn as successive approximations 
to these scalars become available. 

The first of Eqs. ( 1 ) is fundamental to the Gauss method, and standard to any evaluation 
of it is consideration of the analogous equation 

Po = {Zq - C \Z\ - c%Z$)jno, (2) 

where c\ and c 3 are approximated by the ratios of the areas of the corresponding sectors, 

C? = (<3-*0)/(<3-<l) 

c ° 3 = (t 0 -ti)/(h-t 1 ), 


and can be immediately computed. Migration from Eq. (2) to the first of Eqs. (1) requires 
use of the “sector- triangle ratios”, and in view of the well-known power-series expansions for 
these it has been usual to establish this in the form 
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Po = /£>o(l ~7oAo), ( 4 ) 

where ro, the asteroid’s heliocentric distance at time to? i s ®dso related to po geometrically, 
and 70 is another quantity that is modified during the iterative process. 

Although not^essential to the normal orbit- de termination process, it is useful for the 
present purpose also to define a quantity awi 

a° = m 0 /(micf), _ (5) 

and since m 0 and ra a are the sines of the angles between, respectively, the observations at t 0 
and t 3 and at ti and f 3 , in a short-arc or bit the ratio m 0 /mi goes roughly as the time ratio 
c°, in which case q° ~ l.~' - g ~ 1 : ‘ ” ‘7, . - = : ~ 

Although Gauss used Eq. (4) and its geometric counterpart quite rigorously and com- 
puted the sector-triangle ratios in terms of the hypergeometric function, others h„ave instead 
relied heavily on power-series expansions, and since these expan sions ca n certainly divefge, 
Taff’s criticism has been understandable. In the case of a main-belt asteroid observed for up 
to a month or so near opposition, 70 ~ -Rq> where Ro = (Ao + ^ 0 + an ^ convergence 

is to be expected. Nevertheless, acknowledging Taff’s comments, and also because Eq. (4) 
is in reality a pure contrivance, I recommended (Marsden 1985) that the orbit solution be 
carried out by incorporating successive iterates for the ci and c 3 directly into Eqs. (1). 

In the particular situation being considered here, however, the unit vector (/o? mo, no) 
the isolated third observation is unknown. This suggests that Eqs. (2) and (5) should be 
inverted and written 

n 0 = ( Zq - c\Z\ - c%Z 3 )jpo ^ 

rho = amiCj. 

One can therefore assume suitable values of p 0 and a and compute the corresponding n 0 , mo 
and 7 0 = (1 - ml - rag)U The sets of (7 0 ,rh 0 ,no) can then, on the one hand, be transformed 
from the Cunningham system back to the equatorial system and unique sets of values of the 
right ascension and declination, and on the other hand, be substituted for the (/o,m 0 ,no) in 
Eqs. (1) and each time subjected to the full GEM orbit-determination process. The region of 
the sky to be searched can then be restricted on the basis of those orbital solutions that give 
acceptable values of semimajor axes a and eccentricities e (and perhaps also the inclinations 
i, although this requires conversion back from the Cunningham coordinate system). 

It is useful to attempt the computation of a circular orbit from the observations at t\ 
and t 3 . The resulting position at t 0 can be substituted in Eqs. (2) and (5) to give p% and 
a c , say, and the trial values needed in Eqs. (6) are then taken in the vicinity of po = p% 
and a = a c . It is not clear a priori how extended that vicinity should be, although Lnote 
that, for a main-belt object, Eq. (4) implies that po is ~ 0.1 AU larger than po, an d since, 
near opposition, p 0 ~ r 0 - 1 AU, to meet the extremes of perihelion and aphelion distance 
considered by Bowell et al. (1989) requires 0.7 < po ^3.1 AU. 

For an example, I use that of 1985 FZ considered by Bowell et al. (1989), and Fig. 1 
is an adaptation of their Figure. The times of the connected observations are t± = Mar. 
21.301 UT, t 3 = Mar. 24.323, and the isolated plate was taken on t 0 = Apr. 14.287. The 
circular orbit from the March observations yields a = 2.79 AU, p% = 1.924 AU^ c* c = 0.952, 
and Fig. 1 shows that all the acceptable solutions are contained within 0.8 < po < 2.8 AU, 
0.82 < a < 1.02. The original Apr. 14 candidate, denoted by “A”, has an orbit with a = 2.66 
AU, e = 0.206, i = 5°.5. However, Williams (1991) has shown that this observation, but not 
the March observations, refers to a different object. There is another candidate to the north, 
“B”, and this has an orbit with a = 2.36 AU, e = 0.253, i = 4®2. 

If the motion in ecliptic latitude is too high, it may not be possible to calculate a 
circular orbit. Near opposition, and if i is ignored, the radius of the orbit (in AU) can be 
approximated by 


a — 1 + k(t 3 - ti)(l - as)/mi, 


(7) 
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RIGHT ASCENSION 

Fig. 1. The region of the sky (equinox 1950.0) in which 1985 FZ was located on 1985 
Apr. 14, showing the positions in terms of the quantities po and a. The shaded region 
refers to acceptable orbits. The Vaisala locus of perihelic and aphelic orbits (including a 
circle) is shown. “A” is the original candidate, “B” an alternate candidate. 

where the times are measured in days and the object’s mean daily motion s = a~» is in units 
of the earth’s mean daily motion k, i.e., the Gaussian constant 0.017202. For the 1985 FZ 
example, mi = +0.01185, and Eq. (7) can be solved iteratively to give a = 2.73 AU. 

The quantity a describes the foreshortening of the apparent motion, and for a circular 
orbit in the plane of the ecliptic with observations placed symmetrically about opposition 
it follows that 5=1, The foreshortening is greater than or less than unity according as 
to whether the observations are mainly before or mainly after opposition, the approximate 
time of which is given by 

topp = h + Y 3 {t 3 - t 1 )/(Y 1 - Y 3 ). (8) 

For 1985 FZ, Yj = +0.09084, Y s — +0.03992, so that t opp — Mar. 26.692. The foreshortening 
can then be approximated by 

a = (l + i^P) T , (9) 

where, with the time intervals in days, 

T = (t 0 -t 1 )(t opp -t mem ), (10) 

i.e., -86.64 day 2 in the example, f meBn being the mean observation time Mar. 30.304, and 

P — (1 - as 3 )/( 1 - as) + 3a[(l - s)/(a - l)] 2 , (11) 

i.e., 4.12 for a = 2.73 AU, so that 1 + \k 2 P = 1.00061 and a = 0.949. The factor P increases 
between limits of 1 and j as a decreases from oo to 1 AU, or between 3.5 and 5.2 over the 
main belt. Unless |Tj is very large, the expression 1.0006 r is therefore quite sufficient for 
establishing the initial trial value of a for a main-belt object. 
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The new procedure can also be used in the case of Apollo asteroids, for which Vaisala 
orbits are troublesome because these objects can not simultaneously be at opposition and 
perihelion. Circular orbits are also meaningless for Apollo objects, so the relevant ranges of 
p 0 and a must be established by trial and error. The object 1990 MU was observed on two 
nights separated by six days in June 1990. The resulting search region for a night one month 
later extended some 1.5 hr in right ascension, its width increasing from as little as 1° at the 
eastern end to more than 5° at the western. The range of acceptable values of a was ~ 0.3- 
1.3, with po ~ 0.9 AU for small a, up to 0. 8-1.5 AU for moderate a, then down to 0.6-0.9 
AU for large a. Small values of p 0 put the object near the earth and in earth-like orbits, and 
large po involved solutions that made the object distant with large e and i. For the largest 
a the orbits degenerated into Aten type, and the search area became rather nebulous. The 
object was identified at po = 1.28 AU, a = 0.89 with a = 1.62 AU, e = 0.66, i = 24°. 

The case of the Apollo object 1991 JW is interesting because searches were to be made 
on the basis of observations spanning only a single 46- min exposure three weeks earlier. 
Given that the object was retrograding, however, the general character of the solutions was 
similar to that of the 1990 MU solutions, though with slightly larger ranges in po and a not 
anticipated to be greater than 1.1. Although 1991 JW was in the mean time independently 
discovered by another astronomer, the first observer did succeed in recovering it right in a 
corner of the search area. With po — 1.06 AU, dc — 0.45, the object had an exceptionally 
earth-like orbit with a = 1.04 AU, e = 0.11, i = 8°. 

Of course, the procedure outlined in this paper is equally valid however the three ob- 
servations are ordered — if the isolated observation were made during the month before the 
others, it would be reasonable to reverse the order of the t; — and a more detailed write-up 
in the Oct. 1991 issue of the Astronomical Journal discusses how the procedure can be used 
with advantage when a three-observation orbit solution is indeterminate, or where one might 
suspect that an observation is significantly in error. 

The Astronomical Journal write-up also considers how the procedure can be extended 
to handle additional observations. The 1985 FZ example would have fared much better if 
there had been a second night of observations to verify the correct linkage of the object 
in April. Given such a “verification” observation, made at time t Vf r , preferably only a few 
days from t 0 , it is not difficult to compute, still for the time t 0 , a p 0 — p y Q and an a = a v 
that exactly satisfy the verification observation. If the linkage is correct, these should be 
very similar to pg and a°, and the residuals at t 0 and t v „ can then be better distributed by 
adopting the means p$ = %(Po + Po), = H a ° + aV )- Obviously, further observations could 
also be included in these means, and the outcome resembles that of Herget’s (1965) method 
and “poor man’s” least-squares fit, the residuals of the observations at t x and t 3 remaining 
precisely zero. Herget’s example of 1935 QA, with the observations in the order t\ — Aug. 
30, t 0 = Sept. 6, fver = Sept. 23, t 3 = Oct. 21, results in p% = 1.78999 AU, a° = 0.95298, 
pi = 1.78883 AU, a v = 0.95296. The means and a m yield residuals (in arcsec in the 
equatorial system) of +0.4, +0.1 at t 0 and -1.1, -0.6 at t v „. These are a little larger than 
Herget’s values because of the wide separation of to and t ve r . Residuals more comparable to 
Herget’s follow if pg 1 and a m are weighted closer to pj and a v . 
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